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Abstract 

A recently proposed renormalization scheme can be used to deal with nonrelativistic 
potential scattering exhibiting ultraviolet divergence in momentum space. A numerical 
application of this scheme is made in the case of potential scattering with r~ 2 divergence 
for small r, common in molecular and nuclear physics, by the use of cut-offs in momentum 
and configuration spaces. The cut-off is finally removed in terms of a physical observable 
and model-independent result is obtained at low energies. The expected variation of the 
off-shell behavior of the t matrix arising from the renormalization scheme is also discussed. 
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1 Introduction 



Ultraviolet divergences commonly occur in field theories, as a consequence of point-like inter- 
actions in the original Lagrangian, both in perturbative expansions and exact solutions. These 
divergences in perturbative quantum field theory could be eliminated by renormalization to 
yield a scale 0, |[ and with few exceptions the renormalized perturbative series could not be 
summed. 

In quantum scattering and bound states calculations, a similar situation occurs for potentials 
with certain singular behavior at short distances. In momentum space such problems develop 
ultraviolet or large momentum divergence. In these cases, both perturbative and exact methods 
fail to obtain finite solution. 

Motivated by the kind of equations that appear in quantum field theories, the renormal- 
ization ideas have been discussed in configuration || ||] and momentum f| j|, [7|, H spaces for 
potential scattering, where Dirac delta, contact, or zero-range potentials were used. Such po- 
tentials have the advantage to be simple, local and separable, with the additional advantage 
that all discussion can be analytical [|, 0. 

Usually, renormalization is carried out by introducing a sufficiently large (small) cut-off in 
momentum (configuration) space. With this cut-off the divergences of the original problem are 
avoided and a meaningful result obtained. In a successful renormalization scheme the knowledge 
of a physical observable at a particular energy is used to eliminate the cut-off and the final 
renormalized results should be independent of the detailed renormalization scheme in the limit 
when the momentum (configuration) cut-off is removed to infinity (zero). Renormalization of 
these potential models leads to new scales and finite physical observables ||. 

The purpose of this work is to perform numerical renormalization of a quantum mechanical 
problem where analytical solution is not possible. In this way, we can study and understand 
better the on- and off-shell behavior of the resulting t matrix, arising from different regulariza- 
tion schemes. We consider renormalization of S'-wave potential scattering with a local attractive 
potential singular at r — > 0, given by 

V{ r) = V a ?StSl. (!) 

In addition to performing general numerical renormalization with this potential we shall be 
particularly interested in studying the expected off-shell behavior of different regularization 
procedures in the nucleon-nucleon interaction. For a sufficiently large |Vo|, the scattering prob- 
lem with potential (|TJ) exhibits ultraviolet divergence and does not permit solutions ||. At 
negative energies an infinite number of bound states collapse to infinite binding and at positive 
energies the scattering equation has noncompact kernel and does not allow scattering solu- 
tion. For a very small |Vo| these theoretical problems can be avoided but the problem becomes 
difficult to handle numerically. 

It is appropriate to mention that, in the limit of fi = 0, the critical strength given in Ref. [[]] 

is 



2 



where m is the reduced mass of the two-body system separated by distance r. For attractive 
(negative) potentials, with |Vq| > |Vo iCr ft|, one encounters an infinite number of bound states 
with infinite binding. For |Vo| < |Vo )£T it|, these infinite number of bound states are absent. For 
a finite non-zero /x, as the short range divergent behavior of potential (|l|) is unchanged, similar 
properties are expected. This behavior is apparent in our numerical results. 

The motivation behind studying this potential is that similar potentials occur in various 
branches of physics, such as molecular, nuclear and atomic physics f|3T0| | . Also, from a mathe- 



matical point of view, renormalization in quantum mechanics is an interesting field of interest. 
Potential ([!]) has particular interest in molecular and nuclear physics. The meson-exchange 
field-theoretic nucleon-nucleon potentials have terms similar to that in Eq. ([!]). Potentials with 
similar behavior also appear in molecular physics. Such terms are treated phenomenologically 
by introducing an arbitrary cut-off at small r or large momentum so as to cure the original 
divergences and to fit some observables. As the problem is mathematically divergent, the re- 
sults could be very sensitive to the cut-off and it is not clear that if the momentum cut-off 
can be taken to infinity in order to produce a meaningful renormalization of the problem. In 
nucleon-nucleon potentials the finite size of the hadrons sets a cut-off in both momentum and 
configuration spaces, and one should work with a finite cut-off. However, it would be pleasing 
to see that the results do not get drastically modified as the momentum cut-off is taken to 
infinity. 

In the present work the above problem is renormalized numerically by the introduction 
of a momentum or configuration space cut-off in different ways. The cut-off is introduced in 
the Green function or the potential in momentum and configuration spaces, respectively. The 
cut-off is finally eliminated in favor of a physical observable/information, which is chosen to 
be the scattering length and the number of bound states of the system. Eventually, we find 
that by taking the cut-off to infinity, model independent results are obtained. An example 
of this scheme can be found in Ref . ||1 1|| , where the binding energy was used as the physical 
observable. Though in the nucleon-nucleon system one has utmost a single bound state, there 
could be several bound states in atomic/molecular systems. In this study we have used up to 
four bound states in order to explore the applicability of the renormalization scheme to these 
general problems. 

The regularization of the solution by modifying the Green function is standard in particle 
and intermediate energy physics. In atomic, molecular and nuclear physics the traditional 
method is to modify the potential by form factors. This later procedure is comparable to the 
modification of the potential at short distances or large momentum as has been done in this 
work. We demonstrate an equivalence between both schemes for on-shell scattering. 

Though most aspects of the present study is quite general, we shall examine some special 
aspects of interest to nucleon-nucleon scattering. In addition to the on-shell scattering under 
diverse situations, we also studied the problem of off-shell scattering with potential ([!]) with a 
single bound state in some details as that potential simulates nucleon-nucleon scattering. It is 
well-known that the knowledge of all on- and half-on-shell t matrix elements is enough to pro- 
duce the fully off-shell t matrix elements [0. This is why we only study the half-shell function 
|T3|. These elements are also interesting from a physical point of view as they are precisely 
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the elements which can be probed experimentally in a bremsstrahlung experiment [RJ]. From 
the present study we determine the extent of the off-shell variation of the phenomenological 
nucleon-nucleon potentials arising from different regularization schemes. 

In this study, in addition to one bound state, we also considered up to four bound states 
for potential ([!]). The appearance of few bound states is common in atomic and molecular 
physics, though in the nucleon-nucleon system commented above there is at most one bound 
state. The successful renormalization of potential ([!]) in the presence of several bound states 
suggests that the present renormalization scheme should be applicable to similar potentials in 
atomic and molecular physics. 

Similar renormalization has been performed recently in quantum mechanics for the one- 
dimensional x~ 2 potential and the three-dimensional Hulthen potential by Gupta and Rajeev 
[ 15[ , where the interest of such potentials in polymers has been emphasized |L6| . 

The plan of the paper is as follows. In section II, we briefly describe the renormalization 
scheme in quantum mechanics, in case of a general local potential, following Ref. 0. In section 
III we describe the different regularization procedures applied to the scattering equation for 
the potential given by Eq. (|l|). Finally, in section IV, we present results and main conclusions. 



2 Renormalization in Quantum Mechanics 

In operator notation, the Lippmann-Schwinger equation for the t matrix T is given by 

T(E) = V + VG (E)T(E), (3) 

where Gq(E) is the free Green function operator, 

G (E) = (E-H + iO)-\ (4) 

H is the free Hamiltonian and E is the center-of-mass energy. V is the potential that, in the 
present case, is given by Eq. ([!]). 

The partial- wave projection of a matrix element of the operator O (= T or V) in momentum 
space, for a spherically symmetric potential, is given by 

r> oo I 
" 1=0 m=-l 

7T ^ 47T V PI J 

where Y^ m (Q p ) are the usual spherical harmonics normalized to one, Q p represents the two 
polar angles of p, and p = |p|. Pi(x) is the usual Legendre polynomial of order / and Oi(p,q) 
is the coefficient of the expansion. 
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Using Eq.([5|) in Eq.([3]) for a central potential, after the integrations in the angular variables 
(where we use the orthonormalization condition of the spherical harmonics), we obtain the 
corresponding partial wave projection of Eq. @: 

T l (p,p';E) = V l (p,p') + - / q 2 dqV l (p,q)G (q;k 2 )T l (q,p';E), (6) 

7T JO 

where E = k 2 , in units such that h = 2m = 1, with m the reduced mass. The momentum-space 
free Green function is given by 

G (q;k 2 ) = (k 2 -q 2 + z0y 1 . (7) 

With the definition given in Eq.((5|), applying the Fourier transform of the local potential V(r) 
to momentum space, we can easily obtain the momentum space partial wave coefficients, as 
detailed in the following. From Eq.(||), by defining x = p.q/(p?) and applying the usual 
normalization of the Legendre polynomials, we have 

Vi(p,q) = IT 2 jf tfcP,(aO(p|V|q) 

= 7r 2 J ^dxPi(x) ^ J d 3 rV(r) exp (z(p - q).r) 

= ±J\ 2 drV(r) J^dxP^x f^J^ . (8) 



Next, we use the addition theorem given in Eq. (10.1.45) of Ref . [[17 

sin(|p — q|r) 
|p — Ql r 



jo \ r\lP 2 + q 2 - 2pqx = J2( 2n + l )jn(pr)j n (qr)P n (x), (9) 



n=0 



where ji(z) are spherical Bessel functions of order I (see chapter 10 of Ref.|[L7||). Finally, using 
the normalization condition of the Legendre polynomials, we obtain a simplified expression for 
our definition of the /—wave projection of a matrix element of local potential in momentum 
space: 

poo 

Vi(p,q)= r 2 drV(r)ji(pr)ji(qr). (10) 



o 

As we are going to consider only the S'-wave (/ = 0) case, in the following we simplify the 
notation by neglecting the I indices. 

With the above partial wave projection, the S'-wave scattering length a and the off-shell 
function f(p,k) are defined, respectively, by 

a = T(0,0;0), (11) 

and 



In Ref. [0, a constant potential in momentum space, V(p',p) = A, was chosen to exemplify 
the application of the renormalization scheme in Quantum Mechanics. With such potential the 
above equation presents ultraviolet divergence, because the kernel of the integral equation is 
noncompact and Eq. (|6|) does not have scattering solution. 

So, to give meaning to Eq. (|), when it presents ultraviolet divergence, a new regular- 
ized equation should be defined, where the Green function or the potential is replaced by a 
regularized one with a cut-off parameter in momentum space (A) or configuration space (a), 
respectively. Such regularization truncates the large-momentum parts of this equation and 
eliminates the ultraviolet divergence. The cut-off parameter is finally eliminated in terms of 
a physical observable, in order to yield a renormalized t matrix which should be independent 
of the regularization procedure. The Green function and potential regularizations should lead 
to equivalent results, when the cut-off is removed by taking the limit A — > or a — > 0, re- 
spectively. We show in the following that the phase shifts calculated through a Green function 
regularization is very similar to those calculated through a potential regularization in the above 
limit. 

One way of introducing a regularized Green function is to multiply the free Green function 
by a regulator function v(q, A; k) that contains a smooth cut-off A(> k), such that 

The function v(q, A; k) satisfies the conditions 

v(k,A;k) — 1, (14) 

and 

lim v(q,A;k) = 1. (15) 

A^oo 

For an appropriate v(q, A; k) the regularized scattering equation, 

T R (p,p';E) = V(p,p') + - / q 2 dqV{p,q)G R (q,A;k 2 )T R (q,p';E), (16) 

7T JO 



v 2 (q,A;k)= l /A2 7 ; . (17) 



has no ultraviolet divergence. 

In case of a constant potential in momentum space, the following regulator function was 
introduced in Ref. 0: 

(A 2 + k 2 ^ 
(A 2 + q 2 

For other potentials presenting ultraviolet divergence, one can try the same regulator function 
v(q, A; k), or, (in case of stronger divergence) replace it by another convenient function. 

The regularization procedure used above guarantees unitarity because the imaginary part 
of the Green function is unaffected. However, the limit A — > oo reduces the regularized Green 
function to the free Green function and the original ultraviolet divergence reappears. The 
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renormalized t matrix is obtained by eliminating the cut-off A in favor of a physical observable 
at some specific energy 

The above regularization through the Green function, as given in Ref. 0, is on-shell 
equivalent to a regularization through the potential. To demonstrate this, let us regularize Eq. 
(fD using the following regularized potential 



V R (p, p') = v(p, A; k)V(p, p')v(p', A; k). 



;is) 



and keeping the Green function unchanged. 

Now the regularized t-matrix elements satisfy 



2 r°° 

T' R {p,p'-E) = V R {p,p')+- / q 2 dqV R (p,q)G (q;k 2 )T R (q,p';E). 

IT JO 



(19) 



The on-shell equivalence of both schemes can be shown by the relation between T R and T R , 
that can be obtained by replacing V R , in Eq. (|19D , by Eq. (|18|): 



r R (p,p'; E) = v(p, A; k)T R {p,p'- E)v{p\ A; k) 
, v , , r-, <■ , 2 r 2, \ v 2 (q,A;k) 



(20) 



2 f°° 

V{jp,p') + - / q 2 dqV(p,q) 

TC JO 



(k 2 -q 2 + lO) 



T R (q,p';E) 



v{p',A;k). 



Considering that function v(p,A;k) should satisfy Eq. (|14|), all the on-shell scattering 
observables derived by using T R or T' R are identical. Differences between the two approaches 
can only be detected when calculating off-shell observables. 



3 Models for Regularization and Renormalization 

Fredholm reduction of the Lippmann-Schwinger equation for the t matrix has numerical advan- 
tage in dealing with real algebra and one such scheme was used for numerical solution. Auxiliary 
Fredholm real T-matrix equation was solved numerically |Tj| [L9| . The T-matrix equation avoids 
the fixed point singularity and has a much weaker kernel. (For such reduced kernel techniques 



we refer the reader to Refs. |L3| and references therein.) When accurate numerical results 



are needed, the iterative procedure applied to the V matrix integral equation has considerable 



advantages, compared to the matrix inversion methods, as shown in Refs. . 

We apply in this section the renormalization scheme presented in the previous section to 
solve the scattering equation with potential ([!]). For a large enough |Vq|, this potential has 
a singular behavior at the origin and admits a ground state with infinite energy, which is 
physically meaningless. 

To test numerically the renormalization scheme we apply four regularization procedures for 
S'-wave scattering with this potential. In two of them we have introduced a regulator function 
in the Green function as described in the previous section. In the other two approaches we 
introduce a cut-off directly in the configuration space expression of the potential, so that a 
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simple regulator function in momentum space cannot be factorized. We describe the four 
schemes (A, B, C, and D) in the following. 

Choice A: The S'-wave momentum-space representation of potential ([[]), after solving Eq.fllOD 
for I = 0, with help of Eq. (3.947.2) of Ref.f2U|, is given by 



2\xq 2fip fi /i 2 + (p - q) 2 

p arctan — - + q arctan — + - In 



jj, 2 + p 2 — q 2 n 2 + q 2 — p 2 2 jj 2 + (p + q) 2 

with the following limiting behavior 



(21) 



V(0,q) = V(q,0) = —arctan (-) . (22) 

The regulator function in the Green function is taken to be the following cut-off step function: 

v 2 (p,A A ;k) = Q(A A -p), (23) 
where Q(x) = 1 for x > and =0 for x < 0. 

Choice B: The potential V(p,q) is also given by Eq. (21), but in this case the regulator 
function in the Green function is taken to be a smooth cut-off function: 

Ap,A B ;^{f±ff (24) 

Choice C: The regularization procedure is given in the configuration space, through a cut-off 
ac, included in the potential, such that in the limit «c — > the regularized potential reduces 
to Eq. (|^). In this case the potential used in calculation is given by 

r(r + a c ) 

The corresponding S- wave momentum-space matrix element of potential Vc(p, q) is given by 
V e^ ac 7 (z 2 + {p + q) 2 \ 



Vc(p,q) = —z J dzexp(-a c z) In 



2 



Apq J ' \z 2 + (p — q) 

V ln / y 2 + (p + g) 2 \ _ ZVo efiac 1 dz _ z exp(-a c z) 



Aatcpq \H 2 + (p - q) 2 ) ac J [z 2 + (p + q) 2 } [z 2 + (p - q) 2 ] ' 

Vo^o } dy eM-mc/y) ln ( ^ + (p + q )y \ (26) 



Apq J y 2 \/i 2 + (p - q) 2 y 2 J ' 

with the limiting behavior 



V C ( P , 0) = V C (0, p) = / dz <?E^1 = v^ro J dy (27) 



z 2 + p 2 " J jj 2 + (pyY z 
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Choice D: The procedure for regularization is given in configuration space, through a cut-off 
od, such that in the limit — > 0, V reduces to Eq. ([[]). The regularized potential for this 
choice is 

V exp(-fi r) 
(r + a D ) 



V D\ r ) = 72 I I > ( 28 ) 



and the corresponding .S-wave momentum-space potential Vo{p,q) is given by 

z(z — /i) exp(— apz) 



V D (p,q) = 2V e^° Jdz- 



1 



z 2 + {j) + qf}[z^ + {j>-qY} 



2V e" a V J dy ^ + 2 ^ p2 + P J 2) ^ ° p2 V _ q2)2yA - (29) 

For choices A and B the regulator function is introduced in the Green function, as explained 
in section II. For choices C and D, the cut-off is introduced directly in the configuration space 
expression of the potential, before the Fourier transform is taken, such that a regulator function 
cannot be simply factorized in momentum space, as in choices A and B. 

The results for the various renormalization schemes is expected to be independent of the 
detailed regularization procedures. Such expectation is verified in numerical calculation in the 
next section, and its limitation discussed. 



4 Numerical Results and Conclusions 

In the following we present the results for phase-shifts, obtained by using four different reg- 
ularization procedures of the previous section. The reduced mass m has been taken to be 
appropriate for the nucleon-nucleon system; h 2 /(2m) = 41.47 MeV fm 2 . The range parameter 
H of the potential has been taken to be 1 fm" 1 . Though we have chosen nuclear units in our 
calculation, most aspects of this numerical study are of general interest. 

In the renormalization scheme we define the cut-off parameter through one of the choices 
of regularization (for example, choice A), given by Eq. (p3|) , and fix the scattering length to 
the experimental value (a = 5.4 fm, for the S"-wave spin triplet nucleon-nucleon scattering) 
by adjusting the strength Vq of the potential. The scattering length does not fix the number 
of bound states of the system and the strength of the potential fixes this number. In present 
renormalization we considered up to four bound states of the system. For other three regu- 
larization procedures (choices B, C and D) the strength Vq and the corresponding number of 
bound states (obtained with choice A) are kept fixed together with the scattering length by 
redefining the corresponding cut-off parameters. With this procedure we guarantee that we 
are dealing with the renormalization of the same physical interaction via four regularization 
procedures. Presently, the physical observables that we fix by renormalization are the scatter- 
ing length and the number of bound states. The numerical parameters used in calculation are 
given in Table 1. 
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Though we consider up to four bound states of potential (|1|) and also consider the limit of 
very large momentum space cut-off, the situation of interest in the nucleon-nucleon system is 
the case of one bound state, the deuteron, and a finite cut-off determined by the "size" of the 
problem. Hence, for the nucleon-nucleon system the first row in Table 1 is of interest, where 
the momentum space cut-off is set at 10 fm" 1 . We shall explore the off-shell behavior of the 
nucleon-nucleon scattering under this situation. 

In each of Figs. 1-3, we display four sets of curves for phase shift versus energy. In order to 
have compact figures, in all case we set the zero-energy phase shift at it. Each set is indicated 
by an ellipse and labeled by the number of bound states. In each set, there are four curves that 
represent four different regularization procedures. Schemes (A— D) are conveniently indicated 
inside the figures. The cut-off parameters of the four schemes are Aa, Ab, etc, and «d. In the 
renormalization limit the momentum cut-offs and A# should be large, and the configuration 
cut-offs «c and an should be small. Also, the cut-offs should be large compared to the on-shell 
momentum. In schemes A — D this means A\ >> k 2 , A 2 B >> k 2 , a^ 2 >> k 2 , and a^ 2 >> k 2 , 
respectively. 

By comparing the curves in each set of Fig. 1, we observe that the phase shifts are numer- 
ically very close for energies up to 50 MeV, that correspond to k 2 ~ 1 fm~ 2 , or 0.01 A^. The 
deviations start to increase as energies approach the cut-off limit. 

In Fig. 2, we repeat the same plots as in Fig. 1 with larger momentum cut-offs, A^ and 
A B , or reduced configuration cut-offs, ac and ar>. By increasing the momentum cut-off we 
approximate the renormalization limit of large cut-off. In this case, comparing the curves in 
each set, we observe that the phase-shifts are numerically very close for energies up to 500 
MeV, that correspond to k 2 ~ 10 fm~ 2 w 0.025 A^. 

The plots of Figs. 1 and 2 are repeated in Fig. 3 for a further increased momentum cut-off 
(reduced configuration cut-off). As we can see, the phase-shifts, for each set of plots that have 
the same number of bound states, are very close to each other in all extension of the energies 
shown in the graph. The energies can go close to 1000 MeV, or k 2 ~ 20 fm -2 « 0.025A 2 4 . 

As we can see from the three figures, the limit of validity of renormalization increases 
fast with the increase of the momentum cut-off. As expected, the renormalization become 
independent of different regularization procedures, if the momentum cut-off is large enough 
compared to the on shell momentum or equivalently, the configuration cut-off is small enough. 
In this limit the domain of renormalization extends roughly up to a cm. energy E « 0.025-Ea 
where E A = Ti 2 A 2 /{2m). 

For k « fi, pure S- wave prevails and the scattering phase shift decreases linearly as k: 
6 m 7r — ka. This behaviour is consistent with the phase shift of the contact interaction ||: 
tan(<5) = — ka, which for small k reduces to 5 ~ tt — ka. This is clearly represented in Figs. 
1—3. This behavior does not depend on the details of the potentials and on the number of 
bound states. 

From Table 1, and Figs. 1—3, it is interesting to observe that, even for ^ 0, in the case 
of large cut-offs, the critical strength to produce four bound states is consistent with the limit 
given in Eq.(Q). For the nucleon-nucleon system Vq^u = — 10.3675 MeV fm 2 . Note that in 
Fig. 3 we have a cut-off with 30 fm -1 and the corresponding Vq (= —12.482 MeV.fm 2 ). In the 
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limit of very large cut-off and large number of bound states, Vq should approach its critical 
value. 

The cut-off behavior of the renormalization scheme is shown in Fig. 4. In this figure we plot 
the dispersion in phase shifts among the four schemes A — D versus the cut-off parameter in 
the case with four bound states. To show the dispersion in phase shifts we plot the differences 
between the phase shifts AS that show maximum deviation at a fixed energy: choices A and B. 
The three curves correspond to dispersion at three energies 150 MeV, 500 MeV, and 1000 MeV. 
As we can observe in this figure, the dispersion approaches zero as we go to infinite momentum 
cut-off limit, and the results become independent of the regularization procedure that have 
been used. 

We have studied systems that can have up to 4 bound states, but as pointed out before, 
the situation of interest in the nucleon-nucleon system is the first row in Table 1 where there 
is only one bound state and a finite momentum space cut-off of about 10 fm _1 . The situations 
with more than one bound states are of interest in atomic/molecular physics. Of the three 
momentum cut-offs used in this work the one at 10 fm^ 1 seems to be more appropriate for the 
nucleon-nucleon system. 

The present work can be used to conclude about the expected off-shell variation of different 
nucleon-nucleon potential models with the same physical content. Here we have used one 
constraint on the regularization scheme so as to fit the scattering length. In realistic situations 
several constraints are used in phenomenological nucleon-nucleon potentials so that the on-shell 
phase shifts of the different models are practically the same. However, t matrices generated 
from such potentials may have large off-shell variations as pointed out in Refs. |12L |i~4]| . As 



fully off-shell t matrices are determined in terms of the half-shell t matrices |12j we plot in 



Fig. 5 the half-shell function (12) for the various schemes at following energies: E c . m , = 0, 37.5 
MeV, 140 MeV, and 300 MeV. These were the energies explored by Fearing in his study of 
bremsstrahlung |14| with different realistic potentials. 

The half-shell variation in Fig. 5 is not a consequence of the on-shell variations observed 
in Figs. 1—3. At zero energy the four on-shell t matrices are the same but they may have 
distinct off-shell behaviors. In the two cases (choices A and B) where momentum cut-offs 
were employed the regularization schemes are similar and there is no off-shell variation at low 
energies. The same is true for choices C and D. For the parameters used in the calculation we 
note that the momentum space cut-off functions (|23|) and ( 24]) of choices A and B are quite 
different numerically. The same off-shell behavior generated in these two choices suggest that 
such momentum space cut-offs in two phenomenological models should not lead to large off-shell 
variations independent of the values of the parameters used. 

The off-shell variation of Fig. 5 is significantly lower than that observed by Fearing in 
different theoretical models. The larger half-shell variation observed by him was possibly due to 
differences in the physical contents of the phenomenological nucleon-nucleon potential models. 
This indicates different intermediate-range (2 — 4 fm) behaviors of those realistic potentials. 
In the present study the long- and intermediate-range behaviors of the models are kept the 
same, while the short-range part of these potentials are varied in the different regularization 
schemes. This variation should increase as the momentum space cut-off A is reduced. For 
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a fixed A it should increase as the momentum variables k or p of the half-shell function are 
increased. The present half-shell variations noted in Fig. 5 are arising from the inherent 
differences in different regularization schemes. Similar half-shell variation(s) are expected in 
different phenomenological nucleon-nucleon potentials with same physical content. 

In summary, we have tested numerically renormalization schemes using a singular local 
potential, that diverges at r — > as the inverse of r 2 . The scattering integral equation exhibits 
renormalizable ultraviolet divergences in this case. We performed numerical renormalizations 
with four regularization procedures, two with a cut-off in momentum space and two with a 
cut-off in configuration space. We find that the renormalization techniques are applicable to 
this case, such that, for lower energies compared to the cut-off, the results are independent of 
regularization procedures used to regularize the original equation. In the case of the nucleon- 
nucleon system the regularization scheme should maintain a single bound state and employ 
a finite cut-off. In that case we studied the off-shell variation of the different regularization 
schemes and found that the different phenomenological nucleon-nucleon potential models have 
larger off-shell variations than found in the present study. Different intermediate-range 
behaviors of the usual phenomenological potentials are supposed to be responsible for larger 
off-shell variations. The present renormalization in the presence of up to four bound states is 
of relevance in atomic and molecular physics, and suggest that similar potentials in atomic and 
molecular physics can also be successfully renormalized in the present scheme. 

We thank the Conselho Nacional de Desenvolvimento Cientifico e Tecnologico, Fundagao de 
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Table 1: Numerical parameters used in the four regularization schemes. is the regulator 
used in choice A, A# the corresponding regulator of choice B, «c and old are the regulators 
used in the choices C and D respectively. N is the number of bound states. For each set of 
data with the same regulator A^, there is a corresponding figure that is indicated in the first 
column in the left. The scattering length is kept fixed at 5.4 fm in all cases, independent of the 
number of bound states. 



Fig. 


N 


V (MeV fm 2 ) 


A A (fm- 1 ) 


Ab (fm- 1 ) 


a c (fm) 


a D (fm) 




1 


-0.6348 


10 


12.145 


0.020984 


0.008942 


1 


2 


-3.6000 


10 


8.830 


0.131200 


0.044080 




3 


-14.030 


10 


6.488 


0.403600 


0.112400 




4 


-40.132 


10 


5.200 


0.870300 


0.211380 




1 


-0.5306 


20 


25.420 


0.008515 


0.003693 


2 


2 


-2.2325 


20 


19.896 


0.042190 


0.015260 




3 


-7.2221 


20 


15.395 


0.119400 


0.036880 




4 


-17.705 


20 


12.838 


0.238200 


0.065960 




1 


-0.48935 


30 


38.930 


0.005158 


0.002261 


3 


2 


-1.80350 


30 


31.280 


0.022850 


0.008420 




3 


-5.39250 


30 


24.785 


0.062700 


0.020185 




4 


-12.4820 


30 


20.917 


0.122570 


0.035700 



Figure Caption 

Fig. 1 Renormalized phase shifts versus energy using four regularization procedures described 
in the text. The regulator used for choice A is A^ = 10 fm -1 . The corresponding 
regulators for the other three choices, and the corresponding strength Vo that maintains 
the scattering length fixed at 5.4 fm, when varying the number of bound states N, are 
given in the first set of Table 1. 

Fig. 2 The same as in Fig. 1, using the parameters given in the second set of Table 1, that 
have A A = 20 fm -1 . 

Fig. 3 The same as in Fig. 1, using the parameters given in the third set of Table 1, that have 
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A A = 30 far 1 . 



4 The dispersion in phase shifts A5 versus the momentum cut-off A a, for the energies 
150 MeV, 500 MeV, and 1000 MeV, in the case we have four bound states. We define 
the dispersion by the difference between the phase shifts of choices A and B, which are 
enhanced for larger number of bound states N, as shown in Figs. 1-3. 

5 Off-shell function f(p, k) versus p, for different regularization schemes corresponding to 
the first row of Table 1. The curves are labeled by the respective center of mass energies, 
in units of MeV. 
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